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Implementation and Validation of a Wake Model for Vortex-Surface Interactions 
in Low-Speed Forward Flight. 


Background 

•During the crucial phases of vertical take-off, hovering, and low-speed 
forward flight, rotary wing aircraft operate in a flowfield which is signifi- 
cantly affected by the wake of the rotor blades. The predominant feature of 
this wake is the presence of strong vortices shed from the blade tips, which 
then take an approximately helical trajectory, skewed by the motion of the 
craft. The close proximity of these vortices to one another and to the rotor 
blades can significantly modify the performance of the vehicle. To predict the 
performance of rotorcraft, therefore, one must be able to compute the velocity 
field induced by the tip vortices. This is a computationally expensive task, 
since the contributions of every element of the vortex shed from each blade to 
the velocity at each point of interest must be considered, as well as mutual 
and self-induced contributions of the vortices themselves. For a practical 
rotorcraft configuration, with several rotor blades and complex airframes, 
such computations can be prohibitively expensive. 

Beddoes (1) showed that the velocity induced by a helical vortex trajec- 
tory at a point was dominated by contributions due to vortex elements posi- 
tioned at certain "critical" orientations with respect to the point. The 
remainder of the vortex trajectory contributed very little, and the time spent 
in computing these contributions was essentially wasted. Thus he approximated 
each turn of the helical trajectory by two straight line vortex elements 
positioned at critical locations, plus a "remainder" term to account for the 
rest of the spiral. The savings in computer time were dramatic, and Beddoes 
claimed that excellent agreement with observed data could be achieved. 
Beddoes* interest, however, was aimed at the loads on the rotor blade itself 
due to close blade-vortex interactions. 

It was proposed by the author of this report that this method could be 
extended to several useful cases, such as those where the influence of an 
airframe has to be iteratively computed. In addition, other applications have 
since been identified, such as the rapid computation of the induced velocity 
field above and below the rotor, identification and interpretation of 
experimentally observed features, and the validation of more sophisticated 
blade flowfield codes (such as Navier-Stokes formulations) which require the 
coupling of some form of wake model to produce results which can be compared 
with experiment. This work has since been undertaken, and forms the subject of 
this report. 


Objectives 

The objectives for the current year are: 

(1) to code and implement Beddoes* wake model in a VAX workstation environ- 
ment, 

(2) use the code to perform predictions of induced velocity at specified 
points of a rotor wake in forward flight 

and 

(3) compare the results to those obtained using the Scully free wake code. 



As permitted by time, it is also planned to extend the code to enable 
calculation of vortex-induced effects on an airframe model and compare the 
predictions to experimental data. 

Status Summary at the end of six months 

The wake model has been implemented using a VAX 750 and a Microvax II 
workstation. Online graphics capability is available using a DISSPLA graphics 
package. The results shown in the paper by Beddoes have been reproduced 
satisfactorily. The rotor model used by Beddoes has been significantly extend- 
ed to include azimuthal variations due to forward flight and a simplified 
scheme for locating critical points where vortex elements are placed. The 
procedure for computing the remainder term has been modified to remove discon- 
tinuities which became apparent when the results were plotted over the rotor 
disk. The procedure has also been extended to enable computation of velocity 
at points off the rotor blade, and then to enable computation of time-resolved 
and time-averaged velocity at fixed points in the flowfield. The current 
computation procedure has been detailed by Olivier Schreiber in Appendix A, 
and some of the latest results have been shown. 
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A test case has been obtained from current NASA work' ' for validation of 
the predictions of induced velocity. Comparison of results indicates that the 
code requires some more features before satisfactory predictions can be made 
over the whole rotor disk. Specifically, shed vorticity due to the azimuthal 
variation of blade loading must be incorporated into the model. Interactions 
between the vortices shed from the four blades of the model rotor must be 
included. These modifications are being implemented. 

The Scully code for calculating the velocity field is being modified in 
parallel with these efforts to enable comparison with experimental data. 
To-date, some comparisons with flow visualization data obtained at Georgia 
Tech have been performed, and show good agreement for the isolated rotor case. 
Comparison of time-resolved velocity data obtained at Georgia Tech also shows 
good agreement. Modifications are being implemented to enable generation of 
time-averaged results for comparison with NASA data. The efforts to-date on 
the Scully code have been performed under ARO support as part of the research 
on rotor-body interactions, and has not to-date been supported under this 
grant. However, during the next reporting period, the validated Scully code 
will be used for comparison with NASA results and the Beddoes model will be 
incorporated into the Scully code. 

During the next reporting period, efforts will be concentrated on improv- 
ing the modeling of the wake and the rotor. Two test cases are available for 
validation: the NASA data, and the data being acquired at Georgia Tech. Both 
will be used. The code will be revised and updated periodically as changes are 
made and validated. With the campus network reaching completion, it is expect- 
ed to become a routine matter to transmit data and code between Langley and 
Tech. Some such transmissions have been made already. 

Extended Work 

Quantitative flow visualization data recently obtained at Georgia Tech 
now promise to enable modeling and validation of vortex-surface interaction 
phenomena, as discussed in the "extended work" portion of the proposal. This 
will be attempted in the next reporting period. In addition, there is consid- 
erable interest in computing wake effects efficiently to enable accurate 



flowfield predictions in the close vixJriity of the blade using more sophisti- 
cated rotor aerodynamics formulations^ . These avenues are being explored. 
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Computation of induced velocities 

Olivier A . Schreiber 
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This report shows the computer implementation of the calculations of the 
induced velocities produced by a wake model presented in A wake model for 
high, resolution airloads by T.S. Beddoes. 1 ( Reference [l] ) 

Abstract 

An implementation of a trailing vortex system geometry is developed. 

It is combined with a new approach to compute the induced down wash. 

The program development includes added features in addition to the wake 
model presented in reference [l]: Circulation harmonic variation possibil- 
ities, simplified algorithm for the location of the critical elements of the 
wake, computation of downwash outside of the blade and/or outside the 
tip path plane, possibility of including blade bound vorticity induced ve- 
locity with harmonic circulation variation. In addition to the computation 
of induced velocities with wake fixed at a given time, a time averaging 
process over a rotor revolution divided by the number of blades is carried 
out. Thus, results can be compared to experimental field measurements 
obtained in reference [5]. Various options and parameters' influence on 
the results can be investigated in the future and may allow the model to 
give useful results for comparison with experimental data. 

1 Notation 

• ooc Zeroth harmonic coefficient for the bound circulation 

• Q Xc First harmonic, longitudinal, for the bound circulation 

• 6 lf First harmonic, lateral, for the bound circulation 

• b Number of blades 

• c Blade chord 

• c' Blade chord, normalized w.r.t. R 

• c c Vortex core radius, normalized w.r.t. c 

1 Westland Helicopters Ltd. 
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• cp Vortex core radius, normalized w.r.t. R (cp = c c c') 

• Cr Thrust coefficient T/(p7ri? 4 Q 2 ) 

• l biade Blade index ( range: 1 to b ) 

• ix Spiral turn index ( range: 1 to N ) 

• i unit vector parallel to tip path plane and to the vertical plane containing 
the freestream velocity 

• i ac unit tangent vector to the spiral at a given point, [i A c = (dOm/dAi> v )/\\dOm/dAil> v \\) 

• j unit vector parallel to tip path plane and such that j = k x i. 

• 3 AC unit normal vector to the spiral at a given point, (j A c — kx i AC ) 

• k unit vector perpendicular to the tip p?ith plane and directed upward. 

• N Number of spiral turns used 

• PB projection of vector PB on j A c (pb — PB.j A c) 

• r v Radial coordinate of point P, where the induced velocity is being com- 
puted non dimensional w.r.t. R 

• r v Contracted radius, non dimensional w.r.t. R 

• R Rotor radius 

• 7 Z Algebraic curvature radius of projected spiral in tip path plane 

• T Rotor thrust 

• v Induced velocity (positive up) 

• V Free stream velocity 

• y, z Coordinates of a point in the i, j, k vector base. 

• zpB- z ordinate difference between point P and point B. ( zpp — k.PB/R = 
z B - z P — z v - Zp ) 

• a Angle of attack, rotor tip path plane (negative in forward flight) 

• o Radial abscissa integration variable 

• j3 Angle between vortex segment and blade of reference (0 — $ v — ’tfrp) 

• 0 Integration variable 

• A Vv age of a vortex line point. 



• r Circulation strength 

• r Normalized circulation strength (r/(/? 2 Q)) 

• T Corrected circulation strength with respect to vortex core model 

• k Wake skew angle, (- ta.n~ 1 (p x /\)) 

• A Inflow ratio, (p. -f A,), positive up. 

• A, Induced inflow ratio, (v/QR) 

• A,' Reduced induced inflow' ratio, (A./Ajoh ) 

• A t0 Momentum value of induced inflow ratio 

• A l0 h Momentum value of induced inflow ratio at hover ( — \/Ct/ 2) 

• p x Advance ratio, (Kcosa/Oi?) in tip path plane reference frame 

• p, Complementary ratio, (I 7 sin a/QR) in tip path plane, reference frame, 
negative in forward flight or climb. 

• f2 Rotor rotational velocity 

• <3> r Azimuth angle of straight vortex line segment (<£ t , = arg i^c) 

• Vv Azimuth angle of reference blade 

• V’6 Azimuth angle of blade of index i^ude 

• V’P Azimuth angle of point P w'here the induced velocity is computed 

• p Radial distance from vortex center / core radius 

• p Air density, where appropriate 

• <j Rotor solidity, (bc/nR) 


2 Introduction 

A wake model for high resolution airloads (Reference [1]) showed that the prob- 
lem of determining the induced downwash through a rotor at low advance ratios 
can be solved in two successive steps: First, a geometric model of the distorted 
trailing vortex system of the wake is defined and computed from a prescribed 
time averaged downwash distribution. Then, the induced velocity is evaluated 
by positioning simple straight vortex segments at critical locations in the wake 
with respect to the point where the induced velocity is computed and adding 
a remainder term to account for the neglected portions of the wake. The com- 
putation time of the induced velocity at a given point can thus be divided in 
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two parts: The location of the critical positions of the wake (two for each spi- 
ral turn), and the three terms making up the induced velocity for each spiral 
turn (a straight vortex segment for each of the two positions and a remainder 
term). The advantage of this approach is decreased computation time compared 
to a classical approach scheme where each spiral turn is broken down in small 
straight vortex segments which location must also be computed. 

An induced velocity computation is characterized by: One point P where the 
induced downwash is wanted, the rotor position at a certain time as given by 
the azimuth angle of a reference blade Vv - It is possible, then, to compute the 
downwash seen by a blade at a given azimuth if point P is ascribed to lay on the 
blade. It is also possible to compute the induced downwash at a point outside 
the blade and even outside the tip path plane when the rotor is at a particular 
position defined by Vv- However, to compare induced downwash computed 
through numerical models with those measured at experimental facilities, it 
is necessary to perform a time averaging process: The induced downwash at 
a point in the field is computed by placing point P at this point and then 
compute the induced downwash for different positions of the rotor at successive 
increments of azimuth Vv until the rotor has assumed the same orientation in 
the sense of the computation. In effect, the rotor completes one revolution 
divided by the number of blades 6. As a result, the model has been adapted 
for the computation of induced downwash outside the blade and/or outside the 
tip path plane . This, in turn, necessitated simplification of the algorithm to 
compute the location of the critical locations of the trailing vortex system . A 
discontinuity inherent to the former method was thus removed, and with it, 
an extra remainder term, necessary to iron out that spurious effect. A first 
harmonic variation of circulation has been added to represent the variation 
of circulation strength along the tip vortex turns. However, consistency would 
require inclusion of a shed vortex system to assure vortex continuity. A first step 
in this direction has been added in the form of the bound circulation induced 
velocity from the blades. Further study of the model in its present state will 
yield the influence and relative importance of the various options and parameter 
variations available to see wether approximation of experimental results by the 
model can be obtained. 


3 General description of the wake model 

The wake modeled here is the trailing vortex system only. The shed wake is not 
part of the model. The trailing wake is represented by effectively rolled up tip 
vortex lines trailed by each blade from effective vortex radial sources at location 
r Vy contracted radius. Each of these spiral-like lines follow a cylinder skewed 
according to the advance ratio. The model defines the geometry of those lines 
and then computes the induced velocities by approximating each spiral turn by 
a pair of large straight vortex segments positioned at critical points relative to 
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where the induced velocity is wanted. A remainder term for the rest of the spiral 
term is then added. A finite core model is used for the computations as well 
as a radial averaging of the computed downwash for which a realistic sample is 
wanted on a blade. 


4 Description of the wake tip vortex line ge- 
ometry 

A prescribed geometry presented in reference [3] is used. The vertical displace- 
ment of elements of the wake is computed by integration of a prescribed time 
averaged downwash field given in reference [2]. The tip vortex lines are para- 
metrically represented with coordinates x lM r/ v , z v in the tip path plane frame 
of reference moving with the rotor. The parameters are: 

1. ibiade or the index of the blade from which the given vortex line trails ( 
range: 1 to 6 ) 

2. rf*p the angular position of the reference rotor blade (ibiade — 1) ( range: 
0 to one turn ) 

3. A V’,, the age of a vortex line element. It is defined as the angle by which 
the rotor has turned since that element was produced. ( range: 0 to a N 
number of turns ) 


Thus, 


1 

f X v = r t , cos V’t + } 

y v = r v sin rp v l 

(i) 

with 

i. 

U 

II 

' 


= V>r + 27r(i{,!ade - l)/6 - A^t, 

(2) 

It follows from the definition A = fi z + A, that 



- fi z A 4\ + / AjdV’ 

Jo 

(3) 


The averaged downwash field, within the disc defined by the contracted radius 
r v is given by: 

A, = A t o(l + £z'-2?|3/ ,3 |) (4) 

where primed values are normalised by r v and E ~ |k|. To the rear of the disc, 
the averaged downwash field is given by: 

A, = 2A i0 (l - E\y' 3 \) (5) 
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Here, A l0 is the momentum theory value of the induced inflow ratio and is the 
solution of: 

- = [(Vz + A t0 ) 2 + M 2 ] * ( 6 ) 

Three cases present themselves: 

1. cos Vv > 0 The vortex element has always been at the rear of the disc. 


=fi : Arl’ v +2A, # / (1 - £'|y' 3 |)cfV> 

Jo 

z v = » z Arl’ v -+2A, 0 (l-£|y' 3 |)AV>* 


(7) 


2. cos V ; r < 0, *v < ~ r v cos r/) v The vortex element has not moved beyond 
the disc. Thus: 

r&'P* 

z v = /i 2 AV>„ + A, 0 / (1 + Ex' - £;|y' 3 |)<iV’ 

Jo 

=. = M.-AV’r + A,„ / (1 + E(cosi’ v + n x Ail> v /r v ) - E\y' 3 \)drp 

Jo 

z v = fi l Ai/> v + A l0 ( 1 + E(c osip v + n x Aip v /(2r v ) - |y ,3 |))AV’t, (8) 

3. cos < 0, > — t v cos V’* The vortex element has spent part of the time 

within the disc and the remainder to the rear of the disc; thus: 

A V’n 1 A V** 

s. = MzAV, + A, 0 / (1 + Ex' - E\y l3 \)dif’ + 2A,„ / (1 - £|y' 3 |)<ty 

Jo Jk tKi 

where = -2r v cos ip v / fi x thus: z* = fi z Ajp v + A to (l + ^(cosi/^ -f 

y x (-2r v cos V’ t ,/^i)/(2r 1 ,)-|j/' 3 |)(-2r t , cos V’r/Mx)+2A,„(l-.E'|y' 3 |)(AV’„ + 
2r t , cos V’v /Mx) which simplifies to: 

= HzM’v + 2A,„(1 - E\y l3 \)x v / fi x (9) 


5 Computation of induced velocities 

For each spiral turn, the vortex lines trailed by each blade are accounted for, 
in the computation of induced velocity by placing a pair of straight vortex seg- 
ments at critical locations, and adding a remainder term. In this model, the 
vortex lines have constant circulation estimated from the the total thrust. In 
consistency with the simplicity of the model, the computation of the induced 
velocity by the bound vorticity of each blade is easily added as constant circu- 
lation straight vortex segments of length equal to the contracted radius. In case 
the distribution of downwash is wanted along the blade, reference [l] presents 
a radially averaged downwash expression. The radial definition interval used to 
sample the distribution is the length used for averaging. 
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5.1 Pair of vortex segments 

For a given point P where the induced velocity is needed, and for each spiral 
turn of each vortex line, two points B are positioned to support a straight vortex 
segment of length 2 r v and circulation y . The segments are horizontal (i.e. 
oriented parallel to the tip path plane) and are such that their projection on 
the tip path plane is tangent to the projected spiral in the tip path plane. The 
point of tangency is the point for which the normal to the projected spiral passes 
through the projected point of P in the tip path plane. By convention, B is the 
point having the closest projection point to the projection of P in the tip path 
plane. The induced velocity induced by, say segment AC, placed at B, at point 
P can be computed as follows: 


r r ef , r M=c dixMP 

V B — k — — / 

4* J m = a \\M P\\ z 


(10) 


where k is the unit vector perpendicular to the tip path plane and directed 
upward, dl ~ dli A c = Rdai AC , where i AC is the unit tangent vector to the 
projected spiral at point B obtained in the following way: 


IAC 


dOm 

dOm 

dA\p v 


where Om ~ x v i -f y v j. From (1), we obtain: 


dOm 

dArp v 


(in 


Defining j A c unit normal vector to the spiral at point B, directed outward the 
concavity as 

di ac 
d&ip v 

d'AC.. 

we see that for advance ratios smaller than 1, there is no inflexion point on the 
spirals and thus j A c is also: 

3 ac = fc x Me 



T e f f , in algebraic value along the orientation i A c is given in the form based on 
reference [4], 


r eff(p) = r( 



(12) 


F = fo( 1 ~ ai c cos t j> v - b\ c sin il> v ) 


and 


(13) 



where p is the ratio between radial distance and core radius. Reference 1 suggests 
a constant circulation strength given here in non-dimensional form: 


r' — 
1 0 — 


To __ Ct 

R 2Q - b 


(14) 


with a 0 = 2.4 Now, projecting MP on the i A c, j A c, k vector base, 

M P - M B — Rzpsk — RpbJac 

where zp B — k.PB/R — zp — zp = z v — zp , pp ~ P~B.j A c • Then, MP = 
~Rai AC -Rz PB k-Rp B j AC then, dlxM P = dlR{z PB j AC -~PBk) and k(dlxMP) 
-Rpsdl — - R 2 p B da . Also, we have (Af P) 2 = (a 2 + Zp B + p^)P 2 thus: 


r e// r =+r * _ 

4ttP A = _ (a 


Pb doc 

2 + 4b + pD 3/2 


Switching to non-dimensional variables = i>b/(PQ), we 

obtain: 

rii/PB r a=+T ’ _ da 

Ja = -r r (a 2 


iLM 

4;r 


> 2 + 4b+p!) 3/2 


With /3 2 = + p^ and sinh 0 = we obtain 

, r; //PB /■“ g«“h(Tf) 
t,fi ” ' 2it(3 2 J 0 

, _ KuPb 


<w 


cosh* 9 


2tt/ 3 2 tanh ( ar 8 sinh (^)) 


r. 


JfPB 


Thus, we finally have the expression x>^: 

, r ' eJ 1 r v pB 1 


+ Pb) ( z >b + Pb + T v) 1 ^ 2 


with F'ejf = r'(j-^j) and 


(15) 


5.2 Remainder term 


The remainder term accounts for the contribution of the remaining portion of 
the spiral turn. Reference 1 provides the following approximate expression, 
where B represents the closest proximity of the spiral turn to point P: 


= -r: 


0.275 


efI 1 (pm + 4b + ( ! -Pb/t v )zpb)( 1 + 4b) 1/2 


(16) 
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with r' e jf = r'(y^ 7 ) and p 2 = 1 Pn+fru) where p M — 1 -pp/r v if P is outside 

the concavity or pm = 1 if P is inside the concavity (i.e M according to wether 
Pq is negative or positive, respectively.) 


5.3 Blade bound vorticity 

The blades bound vorticity contributions to the induced velocity computed at 
a point P are computed in the form of straight vortex segments of length r tM 
along the blades of the rotor. In the tip path plane frame of reference, with 
origin at the center O of the rotor, the points O, and P have the coordinates 
O : (0,0,0),P : (Rxp, Ryp, Rzp) and for a blade of index hude, the extremity 
C of the straight vortex segment assigned to it has the normalized coordinates 
x c = r v cos^Jm yc = r v sin ^6 and z c - 0 where = rp T 4 2ir(i b i a< ie - 1 )/&. 
The induced velocity induced by segment OC can be computed as follows: 


~ r e// r M=c dlxMP 
4 tt J m ~o || M >|| 3 


(17) 


Now, MP — MO 4 OP — — Raioc 4 R(xpi 4 ypj 4 zpk) and dl ~ Rdocioc 
with i 0 c = cos V’M 4 sin Thus, k(dlx M P) = R 2 (yp cos — xp sin \l> b )da. 

(M P) 2 — {(x P ~ a cos V’b) 2 + {yp - a sin V’&) 2 + tp)R 2 thus: 


r: 


L 


ill 

4r 


(yp cos V’6 — xp sin V’fc) x 
da 


a =o (( x p - acos V'b) 2 + (yp - a sin V’i>) 2 + Zp) 3/2 


(18) 


Now, let A' (a) = (xp — a cos V'i>) 2 + (yp ~ a sin V’t) 2 + z|>- Developping, we have 
A'(a) = a 2 — 2(ip cos 4’b + VP sin i>b) a + ( x p + yj> + zj>) = ABa + Ca 2 for which 
notation, an expression of the above integral can be found as: 



</q _ 2 ( 2 Cr„ + B 

A'(a) - Q ( yj x (r t ,) 


B 

Va 


) 


with Q — 4AC - B 2 Thus, finally, 


T' 


2 2Cr v 4 B 


B 


4 ^( W co S V,-„si„ W ^- 7 _-- a , 


(19) 


with Q — 4 AC - B 2 , C = 1, B = -2(zp cos 4 yp sin ^t>), -4 — x 2 p 4 y|> 4 z 2 p 
and A’ ( r v ) = A 4 Hr t , 4 Cr\ with T' e fj — Y and, 


2 _ {OP x ioc) 2 

1 O'J *> 

R~ c r 
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2 [(a-pt 4 ypj + Zpk) x (cos Vy + sin V’tj)] 2 

' " *’4 

2 _ [(*p s i n V»* - yp cos *F ~p cos ipij — zp sin V»i»i] 2 

p - — — — — 

R~ c r 

2 __ K Jp sin V'h - yp cos V» t ) 2 + -g,] 

P ~ * 2 4 


5.4 Radial averaging of downwash 

The computation of the downwash along the blade involves occurences of close 
blade vortex interactions. This may result in spurious variations of dowmwash 
dependent on radial grid definition and offset. To address this situation, a 
radially averaged expression of the downwash induced by the pairs of straight 
vortex line segments over a radial interval may be used. 

The downwash induced by a straight vortex line segment at point B on point 
P situated on the blade is given from equation 15 as: 


p; P 2 t vPb 1 

1 +p 2 27r (-PB + Pb) ( : pb+Pb + r t 2 ) l/2 


with r' e jj 
Then, 


= and f • Let V = z 7 PB +c% and £ = s* B +r*. 


Vr = 


rVy Pb 

2 tt (p 2 B +v)(p% + £y/2 


( 20 ) 


For an angle (3 between the vortex segment and the blade, the downwash expres- 
sion at the control point P is averaged over the projected length A r* = |Ar.sin/3| 
where A r is the averaging distance, and /? = $*- V’p where <$ r is the azimuth 
angle of the straight vortex line segment, (4> t , — argMc:) Thus: 


r'r,. pgdr 

2vAr ’ y pB _ Ar . /2 (rl + v)( P l + £ yr . 


( 21 ) 


Reference [1] provides the expression resulting from this integral: 

, r'r,. [ (p| 8 +g) 1/a ~* (p|,+g ) 1/2 + * 

" B 2ir2kAr' " [(p^ + £)W + *’ (p^ + f) 1 / 2 - k 


( 22 ) 


where p Bl = p B - Ar'/2, p B , = PB + Ar'/2, A- = (£ - V) 1 ' 7 = (r 2 - cjj) 1 ' 7 , 
and V — Zp B + 4, £ = z% B + r 2 . 
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6 Positionning of the vortex elements 

For each spiral turn trailed by all blades, the position of the two points B have 
to be known by the age of the vortex element they are at so that ps, zpb , 
which are necessary for the computation of induced velocities can be computed. 

The original method involves an empirical formula used for the location of 
the critical elements of the trailing vortex system . However, its implementation 
involves the use of an extra remainder term for the induced velocity computa- 
tion as specified in reference [1] due to a discontinuity in the formulation. Its 
implementation also revealed it to be unsuitable for points where the induced 
velocity is wanted outside the blade and/or outside the tip path plane . The 
method chosen for its simplicity and absolute reliability in all cases is an itera- 
tion method by dichotomy of each spiral turn where the two points B are to be 
found. The criteria for convergence is expressed in terms of the angle between 
the projection on the tip path plane of PB and Jac which, in absolute value is 
ascribed to be smaller than 3 degrees. The algorithm first advances along the 
spiral turn in increments of sixty degrees until the above angle changes sign. 
Then, on the interval thus located, a dichotomy process is performed until the 
tolerance in angle is achieved. The other point B is then found on the remaining 
portion of the spiral turn in the same way. Let a step be defined as a computa- 
tion of the angle between the projections of PB and ;'/c, that is, either after an 
increment along the spiral turn or in a division of the interval where there is a 
change of sign. With the increment of sixty degrees chosen, the number of steps 
necessary for convergence for each point has been minimized to six, in aver- 
age. The advantage of this method is a simplified algorithm, with an adjustable 
tolerance of positionning of the critical elements, at the same amount of com- 
putation time as the original method. No discontinuity effects are brought in 
by the algorithm, and thus the need for the extra remainder term is eliminated. 
However, the main reason for this algorithm is that the original one could not 
handle accurately, points outside the blade, even after extensive modification. 


7 Development of the code 

The development was done in two phases. The first phase saw the experimen- 
tation with the original scheme as presented in reference [l]. The limitations 
of the graphic capabilities of the software available at that time did not allow 
efficient separation of the graphics programming and the computation coding. 
Also, some of the graphics work had to be done on another system where more 
advanced software was available. As soon as adequate software became available 
on the system to be used on a permanent basis, coding of the second version 
of the scheme took place, by taking care of separating the computation core of 
the code from the pre and post processing of the computation. In addition, the 
modularity of the wake model methodology was respected by separating code 
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modules relative to the computation of the critical locations in the trailing vor- 
tex system from those relative to the induced velocity computation. Thus, three 
core modules were created, each compilable, testable independantly and from 
which, validating graphics can be generated: The first module finds the critical 
locations in the tip path plane of the trailing vortex system . The second module 
comprises the subroutine computing the momentum values of induced downwash 
as well as the vertical distorted geometry of the trailing vortex system . Then 
the third module comprises the computation of the induced velocity at a given 
point. The first two modules are tested through one program each, generating 
graphics shown in Figure 1 and 2, respectively. The third module which involves 
computation of induced velocity components is numerically tested through the 
debugger of the compiler. In addition, this last module is used by two ffont-end 
programs producing the numerical results and generating the plots associated 
with those results. The first program computes the induced downwash seen by 
one rotating blade along its span. Figure 3 represents a plot obtained through 
this program. The second program computes the induced downwash field over 
regions of the rotor disc. This program produces both three dimensional and 
contour plots of the results as illustrated by Figure 4 and 5. All programs use 
the same input file and can be linked and run together or at various degrees of 
independance. The variable names fully use the capability of the compiler to 
recognize long names- which are directly analogous to those used in the report 
for good readability of the code. 

8 Results 

The test case treated is the same as the one treated in reference [5]. It is a four 
bladed rotor, at an angle of attack of -3 degrees for an advance ratio fi x = .15. 
From this, a climb ratio p z = Fsin a/QR was computed. The test data is thus 
summarized below: 

• b — 4 

• C'.rp — .0064 

• — .15 

• p z ~ —.0078 

Other input data include the position of the point above the tip path plane 
where the induced downwash is wanted, the vortex core radius, the contracted 
radius and the number of turns used for the trailing vortex system of each blade. 

• z p — .0/7 

• c R = .077 
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• T 


v = 1 . 

• N = 3 

The computation time is of the order of one half second per point where the 
induced velocity by the trailing vortex system is computed. No first harmonic 
variation of circulation was used in the computation of induced downwash . 
Thus, blade bound vorticity induced downwash averages out over time. 

• Figure 1 shows the critical locations in the trailing vortex system as seen 
in projection in the tip path plane for a point P outside the rotor disc. 

• Figure 2 shows the geometry of the distorted trailing vortex system as 
seen in perspective from downstream. 

• Figure 3 shows the plot of the induced downwash seen by a rotating blade 

from azimuth Vv of 0 to 360 degrees along the span from .2 to 1.2 radial po- 
sition (non dimensionalized by R). The results are given in term of A - , re- 
duced induced inflow ratio, (A* = where A, is the induced inflow 

ratio, (v/Q R} a nd A i0 h is the momentum value of induced inflow ratio at 
hover ( — y/Crf 2).' Thus, a positive value corresponds to an induced down- 
wash and a value of 1 corresponds to A’ = A t0 // = — y/Cx/2 = —.056. 
The option of radially averaged induced velocity was used for this case 
producing smoother results than with the simpler expression. 

• Figure 4 shows the plot of the induced downwash field on the region of 
azimuth 0 to 360 degrees along a radial position from .2 to 1.2 (non di- 
mensionalized by R). It can be seen that this plot is very similar to that 
of the induced downwash seen by a rotating blade but it is smoother. 

• Figure 5 shows the contour plot of the induced downwash field on the 
region of azimuth 0 to 360 degrees along a radial position from .2 to 1.2, 
(non dimensionalized by R). This plot contains the same information as 
the preceding plot, induced downwash seen by a rotating blade but it 
is smoother. Overall, there is a good agreement of the magnitude of the 
induced downwash with the experimental results obtained in reference [5]. 
The trends are followed reasonably well on the retreating side and the 
leading edge with a region of upwash extending from 160 to 270 degrees. 
The trends are not followed well at the trailing edge and advancing side. 

9 Conclusion 

The test case computed with the present scheme shows disagreements with ex- 
perimental results obtained in reference [5]. Interesting results can otherwise 
be observed as to the similarity of instantaneous induced downwash seen by a 
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rotating blade and induced downwash field plots. Different options and param- 
eter values can be used to study their influence on the results obtained. Among 
those possibilities are: 

• Utilization of radial averaging and the variation of the averaging interval 
used. 

• Variation of vortex core radius Cr 

• Introduction of first harmonic variation of circulation strength over the 
trailing vortex system and on the blade bound vorticity. 

• Variation of the contracted radius r v 

• Variation of the number of trailing vortex system turns taken into consid- 
eration 

• Variation of the number of increments over which, the induced downwash 
is averaged over azimuth (The present number is 10). 

It is not fully known yet, however what the importance of a shed vortex system 
is and wether it should be added to this scheme. Other factors entering into 
consideration might be lifting surface effects, inboard trailing wake system and 
presence of a body as is the case for the experimental case in reference [5]. 
In conclusion, various options and parameters’ influence on the results can be 
investigated with the present model and may lead to useful interpretations and 
better agreement with experimental data. 
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Figure 1: Critical locations in the trailing vortex system as seen in projec- 
tion in the tip path plane for a point P outside the rotor disc. 









wane view 



Figure 2: Perspective view from downstream of the distorted tip vortex system. 
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Figure 4: Azimuthal distribution of induced downwash in the non-rotating 
reference frame. 




Induced downwash field 



Figure 5: Contour plot of the induced downwash field in the non-rotating 
reference frame. 


